CONSTRUCTION OF iV-BODY TIME-SYMMETRIC INITIAL 
DATA SETS IN GENERAL RELATIVITY 
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Abstract. Given a collection of N asymptotically Euclidean ends with zero 
scalar curvature, we construct a Ricmannian manifold with zero scalar curva- 
ture and one asymptotically Euclidean end, whose boundary has a neighbor- 
hood isometric to the disjoint union of a specified collection of sub-regions of 
the given ends. An application is the construction of time-symmetric solutions 
of the constraint equations which model TV-body initial data. 



1. Introduction. 

The work of Corvino and Schoen [5J lllj shows that one can glue an interior 
region of any asymptotically Euclidean (AE) solution of the vacuum constraints to 
an exterior Kerr or Schwarzschild solution in the following sense: for any chosen 
AE end in any given AE solution of the vacuum Einstein constraint equations, and 
for any chosen interior region (bounded away from infinity) in this solution, there 
is a new AE solution of the vacuum constraints which has an interior region that 
is isometric to the original interior region, and so that the chosen AE end has been 
replaced with an AE end in a space-like slice of a Kerr space-time. If the original 
solution of the constraints is time-symmetric, then the exterior is isometric to an 
AE end in a time-symmetric slice of a Schwarzschild space-time. This work shows, 
remarkably, that one can geometrically shield the details of the gravitational field 
in the interior region, so that sufficiently far away, all that one sees is the effective 
total energy-momentum and total angular momentum of the interior. 

A natural question to pose is whether one can produce solutions of the constraints 
which contain exact copies of the interior regions of two or more solutions. Using 
the gluing techniques of [H [5] [13l [14], one can indeed do this [6j; however the 
solutions resulting from these gluing methods contain two or more asymptotic ends. 
There are gluing results [J] which do incorporate multiple interiors and can produce 
solutions with a single asymptotic end, but all such results to date require that 
the multiple regions be copies of each other, and be placed in certain symmetric 
configurations. 

We show here that in fact one can glue multiple interior regions of solutions of the 
vacuum constraint equations into a solution with a single asymptotic (Schwarzschild) 
end, without imposing any matching conditions on the solutions being glued, and 
without imposing any symmetry restrictions on the placement of the various in- 
terior regions. Our main restriction here is that we consider only time-symmetric 
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solutions of the constraint equations (and therefore Schwarzschild rather than Kerr 
exteriors), and in the solutions constructed, the given interior regions must be suf- 
ficiently far from each other. In future work, we lift the time-symmetric restriction. 
It is not likely that we can weaken the condition that the regions be sufficiently 
distant from each other, or get a more effective bound on the scale of the final 
configuration. 

The motivation for the present work is to set up initial data for the general 
relativistic (Einsteinian) version of the gravitational iV-body problem. We recall 
that one of the signature results of eighteenth century physics is the complete and 
explicit solution of the Newtonian two-body problem for point particles. There 
is no such complete solution for more than two bodies, or for bodies with finite 
extent and interior dynamics, like stars. However, for any number of bodies, and 
for most standard equations of state, choosing iV-body initial data for Newtonian 
gravity involves solving a single (linear) Poisson equation. While the nonlinearity 
and intricacy of the Einstein initial value constraint equations [lj has long suggested 
that choosing iV-body initial data for Einsteinian gravity would be much harder, 
this work (and its followup) indicate that in fact one can set up fairly general 
iV-body initial data in general relativity. 

As noted above, in this paper our focus is on time-symmetric solutions of the 
vacuum constraints, in which case the gravitational field is described completely by 
a Riemannian metric g, which satisfies the single equation R(g) — 0, where R(g) is 
the scalar curvature of g. The aim is to start with a collection of N asymptotically 
Euclidean time-symmetric solutions, and choose a fixed interior region in each, 
thereby specifying each of the "bodies." Then after choosing a collection of points 
which roughly locate the N bodies on a fiducial flat background, one shows that 
there is a solution of the constraints which: i) contains N regions which are isometric 
to the chosen bodies; ii) is a solution of the time-symmetric constraint R(g) = 
everywhere; iii) is identical to a slice of Schwarzschild sufficiently far from the 
bodies; and iv) has the centers of the bodies in a configuration which is a scaled 
version of the chosen configuration. Here the scale factor can be chosen arbitrarily 
above a certain threshold. For example, if N — 2, the distance between the bodies 
can be arbitrarily chosen above a certain threshold value (depending on the ratio 
of the masses of the bodies) . 

To carry out such a construction, we start by recalling that a Riemannian metric 
g on the exterior of a ball in R™ is an asymptotically Euclidean end (to order £) 
provided there are coordinates in which, for multi-indices < £, 

(i.i) -*«)(*)! = o(\x\-M-W) , 

where d is a partial derivative operator. We require £ > 3 in what follows, and we 
take g to be smooth, at least C 4,a for some a £ (0,1). A complete Riemannian 
manifold {M, g) is asymptotically Euclidean if there is a compact set K C M so 
that M \ K is the disjoint union of (finitely many) AE ends. 
We now state the main theorem. 

Theorem 1.1. For each k = 1,...,N, let (Ek,g k ) be an n-dimensional (n > 3 
AE end with zero scalar curvature and ADM mass nik > 0, and let Uk C E^ be 

N 

a pre-compact neighborhood of the boundary dEk- Let M cxt = M™ \ (J B^, where 

fe=i 

£?i, . . . , B 7v are pairwise disjoint closed balls in W l . Then for each e > 0, there is a 
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Riemannian metric g t on M GX t of zero scalar curvature, so that (M ex t,g e ) contains 
a neighborhood of the boundary dM cy ± which is isometric to the disjoint union 

N 

U (Uk,g k ), and (Afextjffe) has one AE end which is identical to a Schwarzschild 

k=l 

k 



exterior with mass m{g e ) satisfying m(g e ) — ^ m k 



< e. 



The theorem allows us to glue together any finite number of AE ends which solve 
the time-symmetric vacuum constraint equation (i.e., which have vanishing scalar 
curvature) into a solution with one AE end. Moreover, the construction is local 
near infinity in each end E k , i.e., any given compact subset of E k can be realized 
isometrically in the final metric (M,g e ). As a direct consequence, then, we can 
take a collection of N Riemannian n- manifolds (M k ,g k ) in each of which we have 
chosen an AE end Ek of zero scalar curvature and positive ADM mass m k , we let 

N 

M ext = R™ \ (J Bk as above, and we can produce a Riemannian manifold (M,g e ) 
fe=i 

with the following properties: M is diffeomorphic to the quotient of the disjoint 

union M ext U ( JJ M k \ E k \ where we identify dE k and dB k ; g e = g k on M k \ E k ; 

(Mext,p e ) has vanishing scalar curvature and one AE end, with ADM mass m(g e ) 
k 



satisfying 



< e. The result of this kind of construction may be 



i=l 

interpreted as an iV-body initial data set. There is some flexibility from the fact 
that the construction only deforms the original solutions near infinity in the chosen 
ends. For instance, we can allow the N original solutions (bodies) to have multiple 
ends, and we can allow nonzero scalar curvature supported away from the chosen 
ends, which might be interpreted as energy density for matter or a field inside the 
bodies. As one of the goals is to glue multiple AE ends into a single AE end, we 
state the following corollary, which follows directly from Theorem ll.il note that if 
each of the M k below is diffeomorphic to R n , then from the above discussion, we 
can choose M to be diffeomorphic to R™ as well. 

Corollary 1.1. Let (M k ,g k ), n > 3, k — 1,...,N, be complete AE solutions 
of the time- symmetric constraint R(g k ) = 0, each of which has one AE end U.lJjj , 
with ADM mass m k > 0. Let U k C M k be chosen compactly contained sub-domains. 
Then for each e > 0, there is a complete Riemannian manifold (M,g e ) with R(g e ) = 

N 

0, which contains a region isometric to the disjoint union (J (U k ,g k ), and with only 

fe=i 

one AE end which is identical to a Schwarzschild end with mass m(g e ) satisfying 
k 

TO (5e) ~ J2 m k < e- 



2. Preliminaries 

We begin with a proposition which follows from [3] or [TU]. The metrics intro- 
duced here give a family of metrics into which the ends will be glued to prove the 
main theorem. 

Proposition 2.1. There is a family of metrics 7t, t € (0,<5), on R" of zero scalar 
curvature which converge as t — > + to the Euclidean metric in C k (for any k ), and 
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are exactly Schwarzschild on {x : \x\ > 1} of positive mass m(7t) and centered at 
x = 0, with lim mOjt) = 0. 

We also recall two results from [8] , cf. [11] . Let K be the span of the set of 
linear and constant functions K := span{l, x , . . . , x n }. 

Proposition 2.2. Let 51 C R" be a bounded smooth domain, and let C S C^°(^) 
be a cutoff function. There is a neighborhood U of the flat metric in C 4,a (f2) so 
that for any D,q CC 0, there is an eo > so that for smooth metrics g € U and 
functions S € C^°(f2o) with \\S\\c°.<* < eo, there is a smooth metric g + h so that h 
is supported in fl and satisfies a bound \\h\\ C 2, a < C\\S\\c°.<*> and with 

R(g + h)- (R(g) + S) e (K. 

The map (g, S) <—> h is continuous. 

We will see in the proof of Theorem 1 1.1 1 in the next section how this proposition 
is a key ingredient in the proof of the following theorem. 

Theorem 2.1. Let [M,g) be an n-manifold with zero scalar curvature containing 
an AE end £ of ADM mass m(g) ^= 0, and let £ r C £ be the region corresponding 
to {x : \x\ > r} in AE coordinates. Let k be a nonnegative integer. Then for any 
e > 0, there is an R > and a (smooth) metric g with zero scalar curvature and 
11$ — 9\\c k (£) < e (norm measured with respect to the Euclidean metric in the AE 
coordinate chart), with \m{g) — m(g)\ < e, and so that g is equal to g on M \ £r, 
and g is identical to an AE end of a standard Schwarzschild slice on £^r ■ 

Remark 2.2. This statement is a slight generalization of the result in [S], in which 

g had the following form on £: gij(x) = (l + ( n _i™ x \ n -i ) 4 ^" 2 ^Sij + hij(x), with 

\df^hij(x)\ = 0(|:r| - '' 3 '~(™~ 1 )) for \(3\ < 3. Under this condition, or more generally 
assuming that g has enough approximate parity symmetry so that the center of 
mass is well-defined [51 [TT] , then the center of the Schwarzschild end can be chosen 
close to that of {£,g). The proof of Theorem 12.11 is a trivial modification of the 
original argument. 

3. Proof of the main theorem 

We now use the results assembled in the preceding section to prove the main 
theorem. We will prove the result in dimension n = 3; the argument in the other 
dimensions n > 3 is essentially identical. 

N 

Proof of Theorem \l.l\ Choose c\, . . . , cjv € K 3 , with ^2 m kCk — and with \c^\ > 

k=\ 

5. Moreover we arrange the closed balls {x : \x — Cfc| < 4} to be pairwise disjoint. 
Let Bo C t 3 be a Euclidean ball of radius Rq = 5 + max{|ci|, . . . , |c/v|}, centered 

N 

at the origin. Let mj- = ^ m^. 

fe=l 

( N \ 

We now construct a family of metrics 5( £iAliC ) on Bq \ I 1J {x : \x — c^\ < 1} ) , 

for < e <C 1, 1/2 < /i < 3/2, and eel 3 with \c\ < 1. Let j e = j te be one of the 
metrics from the Proposition 12. 1[ with mass m{^ t ) = emr, and let tp be a smooth 
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nondecreasing function so that ip(t) = for t < 9/4 and ip(t) = 1 for t > 11/4. 
Define g( e , M , c ) on B \ y {J {x : \x — c k \ < 1}J as follows: 



.4 



• on each = {x : 1 < \x - c k \ < 2}, g {e ^ <c) (x) = (1 + gj^j ) SEuci(x) 

• on each = {x : 2 < \x — c k \ < 3}, 

9{e,^c){x) = (l-^(|a;-c fc |))(l+2|f^ r |) 4 5Euci(x)+V'(|a;-c fc |)7 £ (Ai- 1 (a;-c)) 



AT 

• on Bq 



h\( \J {x : \x ~ c k \ < 3}Y <?( £>Al ,c)(aO = 7«(m ~ c ))- 
v fe=l 7 

/ N \ 
Let O = Bq \ I IJ {x : |x — Cfc| < 2} J . Let (bea smooth cutoff function, which 

^ fe=i ' 

we can define as follows: 

, N , 

• C(x) = 1 on {x : \x\ < Rq - 1} \ U {% ■ \x - c k \ < 3} 

v fe=i ' 

• = V'G 1 ~ c k\) on A k = {x : 2 < \x ~ c k \ < 3} 

• C(x) = ip{R + 2- \x\) on {x : R - 1 < |x| < i? }- 

There is an e\ > so that for all (e, ^i, c) e 6 = {(e,/i,c) e 1 x 1 x I 3 : < 
e < ei, 1/2 < < 3/2, |c| < 1}, g^^^c) is in the neighborhood U from Proposition 
E2 ||3(e,/i,c) - .gEuci||c 5 < Ce, and moreover \\R{g( t ^.c))\\c°.» < eo- Thus there is a 
constant C > so that for any 8 := (e, /i, c) £ 0, there is a small deformation ft, 
(obtained by applying Proposition 12 . 21 with g — gu <ll<c ) and S = —R(g( t ,fj,,c))) with 
the following properties: h e is smooth on Bq and supported on f2, R{g( e4l , c ) + h e ) = 

3 3 

a,j(xJ , where x° := 1, and ||/i e ||c 3 (n) + S l a jl — Furthermore, given e, each 

dj is continuous in (fi,c). We argue that for e small enough, by choosing fj, and c 
appropriately, we can arrange for the integrals J fi s'i?(9( ei(1 . c ) + h e ) dx to vanish 
for I = 0, 1, 2, 3. From this it follows immediately that a,j = for j = 0, 1, 2, 3, and 
hence R(g( e ^, c ) + h e ) — 0, as desired. 

Let g( e ^, c ) — g 8 f° r notational convenience, and let h = g e — <?Euci- Note that 
we can choose C > so that ||/i e ||c 3 (n) < Ce for all 9 £ Q. 

3 

NowletL(Zi) = —A(tr(h))+div(div(h)) = (hij,ij—hu,jj) be the linearization 

at the Euclidean metric of the scalar curvature operator g i— > i?(<?). Note that 
L*(x £ ) = for £ = 0,1,2,3. We also note the coordinate formula for the scalar 
curvature 

3 

3 

(3.1) = E 9 lJ g km {g t k, jm -g l3 , km ) + Q(dg), 

where Q(dg) is a quadratic function in the partials dg contracted with the metric 
g. We can apply (|3.ip to the metric g — g e + h e = <?Eucl + h e + h e to obtain the 
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following: there is a C\ > so that for small e, we have 



(3-2) = E - 9ii,jj) + L(h e ) + B(h e , h e ) + Q(h e ), 



where 



< C 1 \\h^ C '<n)-\\h^c><n)<C 1 C a <?, 



|Q(^)| < Cill^H^.^ < Cl C 2 e 2 



c 2 {n) ■ II" llc 2 (n) 

2 

C 2 (f2) 



We also note that for a Schwarzschild metric g s of mass m and center c, i.e. 

3 



1+ 2l^y ) 5Euci, the leading order term in the expansion £ (5?.^ - 5^) 1/ 
(where i/ e = t^t is the Euclidean unit normal to the sphere |x| = r) is just 



3 / \ 3 



3 \ -rJ 



1 m \ 1 m \ ^ [ x c x J \ a;- 

- 4 V 1 + W\J \W) ^ \W\ ij " W R" 

Let c?/i e be the Euclidean surface measure. Then there arc functions Fe(m, c,r), 
with |i^(m, c, r)| < C2m 2 /r, where C2 is independent of |c| < r/2, so that 

(3.3) y E ~9u,j) v l dfi e = 16ttto + F Q (m,c,r), 



and for £=1,2, 3, 



3 3 
(3.4) 



lQirmc 1 + Fg(m, c, r). 



Using integration by parts along with (13. 2\i and (|3.3p , and the fact that h and 
its partial derivatives vanish at dfl, we obtain 



/ R{g e + h 6 ) dx = 167re(^ - l)m T + 0(e 2 
Jo 



Let Efc = {x : \x — ct\ = 2} and S = dB = d{\x\ < Ro}, and let Co = 0. 
We now use (|3.2[) , integration by parts along with L*(x e ) — 0, as well as (|3.3I) and 
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(gill) , to obtain for £=1,2, 3, 

3 

/ x*R(g 6 + h e ) dx = I x l J2 (&« " + °( £2 ) 

Jn Jn „• „•_-, 



N 3 3 

E 



+167r^em fe 4 + 0(e 2 
fe=l 

3 



/ (** E (4* - sLjH - E(~9>i - tit))** + o{^) 

■J A A — 1 „- — 1 



— l i—X 

16Tre^m T c e + 0(e 2 



N 

where we have used the condition m kCk — 0. 

fe=i 

Thus we see by continuity that for e small enough, for some /i and c, with 
(/i — 1) = 0(e) and c = 0(e), that we can arrange the above integrals to vanish. 
This means for 9 — (e,/i, c), the resulting metric g e := + h? is scalar-flat and 
is exactly Schwarzschild (with the respective masses) near dfl, so we extend it to 
R 3 \ Bq using the exterior Schwarzschild metric. 

To finish the proof of the theorem, we note that by applying Theorem 12.11 we 
can assume without loss of generality that for each fc, there is an Rk > 1 and there 
is a coordinate system x for Ek, so that on {|x| > Rk}, g k has the Schwarzschild 
form g k (x) = (l + <Jy. Indeed, in applying the construction in [5] to the given 
metrics in Theorem ! 1.1[ the original masses can be made to change by an arbitrarily 
small amount (in our case by less than say), by taking Rk = 2R sufficiently 
large; we can assume as well that |x| < R on each Uk- We assume this has been 
done and the masses have been re-labelled to mi, . . . , mjy- 

, N \ 

Let O a = R 3 \ I (J {x : \x — Ck/a\ < — } )■ We can now scale the metric g t by 

pulling back a scaling of coordinates f e : Q e — > fii, f e (%) = ex, for e < In the 
resulting metric g e = ^r/ £ *<?e, the masses of the Schwarzschild neighborhoods of the 
boundaries \x — Cfc/e| = - > 2R are mi, TO2, ■ ■ ■ , mjv, and near the outer boundary 
is firriT- Given eo > 0, we can choose e small enough so that |m(g e ) — nvr\ — 
— 1)tot| < eo- Note that the resulting configuration of centers is the original 
configuration c\ , . . . , cjy scaled by j . Finally we can glue in the metrics g k on 
(E k n{x : \x\ < f}) D U k isometrically, by identifying {x : \ < \x-c k /e\ < |} C £l e 
with E k n {x : \ < \x\ < f }. □ 

Remark 3.1. In the above construction of g e and g e = g 9 + h e , we could allow 
one or more of the m,k to be non-positive. The proof above uses the positive 
mass metrics from Proposition 12. 1[ which were used to interpolate the various 
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Schwarzschild metrics about each Ck, so this requires tot > 0. One could instead 
interpolate using the Euclidean metric if one also adds another gluing annulus near 
the boundary dP>o to patch the Euclidean metric to an exterior Schwarzschild, and 
the proof follows mutatis mutandis. Using this variant of the construction, we could 
in principle allow cases with tot < 0. 



A natural question to ask is whether the TV-body configurations we produce con- 
tain apparent horizons (minimal spheres), apart from any minimal spheres already 
present in the bodies. We want to know whether we introduce any new horizons in 
the construction, which would shield one or more of the N bodies from the others; 
we also note that the construction in Theorem 12.11 is a small perturbation near 
infinity in an AE end, and so it does not introduce any new minimal surfaces into 
the geometry. In this section we show that for n = 3 and for small enough e, no 
minimal surfaces intersect This can be proved by various means (cf. [7]), and 
one such way is to use the following proposition from [9] , which relies on the Rie- 
mannian Penrose Inequality [2] [12] , cf. [3] . For completeness we record the proof 
here. 

Proposition 4.1. Let (M,g) be a complete AE three-manifold with nonnegative 
scalar curvature. Suppose that the sectional curvatures are bounded above by C > 0. 
Then if the ADM mass m satisfies m\fC < ^, there are no closed minimal surfaces 



Proof. If there were a minimal surface in M, then by the analysis in [12] (which uses 
the convex barrier spheres near infinity in an end) , there would exist an outermost 
minimal sphere E. Let A(E) denote its area. 

We note that the Gauss equation implies an upper bound on the curvature of the 
horizon E as follows. Let p be a point in E, and let {ei, e 2 } be a basis of Tp(E) in 
which the second fundamental form II is diagonal. Let k\ =II(ei, ei), k 2 =II(e 2 , e 2 ) 
be the principal curvatures of E at p. Since E is minimal, we have k,\ + k 2 = 0, and 
so K\tt2 < 0. The Gauss equation then yields 



where R denotes the curvature tensor of £, and R the curvature tensor of M. 

We now invoke the Gauss-Bonnet theorem, which together with the preceding 
inequality yields 



4. Apparent Horizons 



in M. 



= R(ei,e 2 , ei, e 2 ) = R{e\, e 2 , e\, e 2 ) + kik 2 < C 




Combining this with the Penrose Inequality to > 



A(S) 



we obtain 




□ 



We now consider complete AE solutions (Mk,g k ), n > 3, k = 1, . . . , N, of the 
time-symmetric constraint R(g k ) = 0, in each of which we have chosen an AE end 
Ek and compactly contained sub-domains Uk C E^. We apply the construction in 
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the proof of Theorem ll.il to obtain a family of AE metrics (M,g e ); we can write, 
using above notation, M = £l' e U V e , where 

N 



Recall Uk c {x G B^. : < i?}, and we have let Ty. be the transition region {x £ 
Ef. : R < \x\ < -), comprised of the annulus where the original metric transitions 
to the appropriate Schwarzschild metric (as in Theorem l2.ip . together with another 
annulus on which the metric is Schwarzschild. We note that g e converges as e — > + 
to the Euclidean metric on f2 e , and each transition region is foliated by convex 
spheres. 

Proposition 4.2. For sufficiently small e, any closed minimal surface in (M,g e ) 

N 

must be contained inV e \ [J Tk- 

k=l 

Proof. If there were a minimal surface whose intersection with M\V e were nonempty, 
then applying arguments in [12] , there is an outermost minimal sphere E £ with re- 
spect to the end contained in M ext into which we have glued V e , with E n (M \ V e ) 

nonempty. Let = V e \ I U {x £ E}. : i < \x — Ck\ < f } ) ■ We re-scale the met- 



ric by a constant so that on f2 £ the metric is /*g e , which is uniformly close to the 
Euclidean metric. Now it follows by applying the argument in Proposition ^. II that 
for small e there cannot be an outermost horizon contained in n e . Thus E e will 
clearly have to intersect both V e ' and M \V e . By the Penrose inequality, and the 
estimate of /i, there is a constant C2 so that A(Y^ e ) < C2e 2 . 

On the other hand, let p £ E e fl (M \ V £ ). There is an r £ (0,1) and for 
each e an r e £ (ro, 1) so that if B^ 1 (p) is the geodesic ball of radius r about 
p, then S = E e n B^{p) is a smooth surface with nonempty boundary (since 
Y> e C\V^ is nonempty). We can now apply Schoen's curvature estimates for the stable 
hypersurface S [15]: there is a C so that for all small e, the second fundamental form 
on E e n B^/ 2 (p) has norm bounded by C. Hence near p, E e n B^^ 2 (p) is graphical 
with bounded gradient and curvature over a ball of fixed size in T p E e . The area 
of E e is thus uniformly bounded from below, either by the graphical description 
or volume comparison (since the intrinsic ball Bf^ 2 (p) C E e n B^^ 2 {p)). This 
contradicts the area estimate from the Penrose Inequality. 

Thus E £ C V e . But then by the maximum principle, it follows that E e C V £ , 
since dQ' e n V e is a union of N convex spheres. Similarly, E e cannot intersect any 
of the transition regions. □ 
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